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Abstract 

In this article, we assume the Z c ( 4200) as the color octet-octet type axial-vector molecule¬ 
like state, and construct the color octet-octet type axial-vector current to study its mass and 
width with the QCD sum rules. The numerical values M Zc ( 4200 ) = 4-19 ± 0.08 GeV and 
Fz c ( 42 oo) w 334 MeV are consistent with the experimental data M Zc ( 4 2 oo) = 4196t2sAi3 MeV 
and r Zc ( 42 oo) = 370170^^32 MeV, and support assigning the Z c (4200) to be the color octet- 
octet type molecule-like state with J pc = l -1 . Furthermore, we discuss the possible as¬ 

signments of the Z c ( 3900), Z c ( 4200) and Z(4430) as the diquark-antidiquark type tetraquark 
states with J PC = l -1 . 

PACS number: 12.39.Mk, 12.38.Lg 
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1 Introduction 

In 2014, the Belle collaboration analyzed the B° -A K~n + J/il> decays with the full T(45 l ) data 
sample corresponds to 711fb~ 1 data sample collected by the Belle detector at the asymmetric- 
energy e + e _ collider, and observed a resonance (named Z c ( 4200)) in the mass spectrum 

with a statistical significance of more than 6.2 < 7 , the measured mass and width are M Za ^ 4200 ) = 
41961^1^ MeV and T Zc ( 4200 ) = 370tyol™2 MeV, respectively [US]- The preferred assignment 
of the quantum numbers is J p = 1"*". 

In 2007, the Belle collaboration observed a distinct peak in the mass spectrum in the 

B° —> decays with the statistical significance of 6.5 a [3|. In 2014, the LHCb collaboration 

analyzed the B° —> K + decays by performing a four-dimensional fit of the decay amplitude, 

and provided the first independent confirmation of the Z( 4430) _ and established its spin-parity to 
be J p = 1 + dj. 

In 2013, the BESIII collaboration studied the process e + e _ —> and observed a 

structure Z c (3900) in the mass spectrum [5], Then the Belle and CLEO collaborations 

confirmed the existence of the Z c (3900) [SEI- Although the quantum numbers are not measured, 
the assignment J p = 1 + is favored if the decay (3900) =*= -a takes place in relative iS-wave. 

In 2014, the BESIII collaboration studied the process e + e _ —> nDD* and observed a distinct 
charged structure Z c ( 3885) in the (DD*^ mass spectrum, the assignment J p = 1 + is favored 
[8]. We tentatively identify the Z c ( 3900) and Z c (3885) as the same particle according to the 
uncertainties of the masses and widths. For more literatures on the X , Y, Z mesons, one can 
consult the recent review [5], 

The quark constituents of the Z c (3900), Z c (4200) and if(4430) are ccud or ccdu , there are three 
analogous decays, 


Z c (3900) ± -A J/V’7T ± , 

Z c (4200) ± -A J/i/jtt± , 

if (4430) ± -A V'A ± , (1) 

which take place through fall-apart mechanism. The mass differences are Af^( 443 o) ~~ Mz c ( 3900 ) = 

576 MeV and — Mj= 589 MeV, so it is natural to assign the Z{ 4430) to be the first radial 
excitation of the Z c (3900) (TDj [TTJ T2] , Naively, we expect the tetraquark states have large decay 
widths, the Z c (4200) and Z( 4430) are good candidates of the tetraquark states according to the 
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widths r z c ( 420 o) = 370i7ol™ 2 MeV [Tj and T Z (aazo) = (172 ± 13 + 34 ) MeV g]. In Ref.[I3], we 
study the axial-vector hidden charm (and hidden bottom) tetraquark states in details with the 
QCD sum rules and obtain the value M Z ( 4430 ) = (4.44 ± 0.19) GeV. In Ref.[Hj, Chen and Zhu 
study the vector and axial-vector charmonium-like tetraquark states with the QCD sum rules and 
obtain the value Mjp =1 + = (4.16 ± 0.10) GeV. In Ref.[15], Chen et al assume the Z c ( 4200) to 
be the axial-vector tetraquark state and calculate its decay width with the QCD sum rules. In 
Ref. .13), the QCD spectral densities are calculated at the special energy scale fj, = 1 GeV, while in 
Ref.[14j, the QCD spectral densities are calculated by taking the MS masses m c {m c ) and choosing 
the vacuum condensates at the energy scale /r = 1 GeV. If the Z c ( 3900), Z c ( 4200) and Z( 4430) 
have the same quantum numbers J PC = l -1 , they cannot all be the ground state axial-vector 

tetraquark state with the same quark structure. 

The X , Y and Z mesons have been studied extensively by the QCD sum rules HMMSEH, 
but the energy scale dependence of the QCD spectral densities are not studied. In previous works 
dans hhi mi [ 201 , we explore the energy scale dependence of the QCD sum rules for the hidden 
charmed and hidden bottom tetraquark states and molecular states in details for the first time, 
and suggest a formula, 


d = yjMx/r/z - (2Mq) 2 , (2) 

with the effective masses Mq to determine the energy scales of the QCD spectral densities. 

The quarks have color SU( 3) symmetry, we can construct the tetraquark states according to 
the routine quark —> diqiurk —> tetraquark, 

(3® 3) ® (3® 3) = (3® 6) ® (3® 6) = 3® 3® ••• = 1 © 8® ••• , (3) 


or construct the molecular states and molecule-like states according to the routine quark —> 
meson (or meson — like state) molecular state (or molecule — like state), 

(3® 3)® (3® 3) = (1 ©8) ® (1 ©8) = (1 ® 1) © (8® 8) © ••• = 1 © 1© •• • , (4) 


where the 1, 3 (3), 6 (6) and 8 denote the color singlet, triplet (antitriplet), sextet (antisextet) and 
octet, respectively. 

In the scenario of tetraquark states, we study the 3 ® 3 type (or the diquark-antidiquark 
type) scalar, vector, axial-vector, tensor hidden charmed tetraquark states and axial-vector hidden 
bottom tetraquark states systematically with the QCD sum rules [12, 17 1 , [18], and assign the 
V(3872), Z c (3900/3885), Z c (4020/4025), F(4140), Z( 4430), Y(4660/4630) and .2/(10610/10650) 
to be tetraquark states tentatively, 


V(3872) = 

i ([ cu]a[cu}s + [cd] A [cd]s + | 

cu\s[cu] A + [a 

Z c (3900/3885) = 

—7= ([ cu] A [cd]s ~ [cu\ s [cd\ A ) 

(with l -1 ), 

Z c (4020/4025) = 
Y (4140) = 

[cu]u[cc?]a (with 1 + or 2 ++ ), 

[cs]a[cs]u (with 2 ++ ), 

Z(4430) = 

-j= ([cu] A [cd\ s - [cu\s[cd\ A ) 

(with l -1 ), 

F(4660/4630) = 

-j= {[cs] A [cs]p - [cs] p [cs]a) 

(with 1 ), 

Z b ( 10610) = 

-i= ([bu] A [bd\ s - HsHa) 

(with 1^ ) , 

.2/(10650) = 

[bu] A [bd] A (with l +_ ), 



(with 1 ++ ), 
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( 5 ) 



where the [Q<z]s, [Qq\A and [Qq]p denote the scalar, axial-vector and pseudoscalar diquarks in the 
color antitriplet 3, q = u : d,s and Q = c,b. 

In the scenario of molecular states, we study the meson-meson type (or the 1 g> 1 type) scalar, 
axial-vector and tensor hadronic molecular states with the QCD sum rules [ESSO], and assign the 
0(3872), 0(3900/3885), V(4140), 0(4020/4025) and 0(10610/10650) to be the molecular states 
tentatively, 


X(3872) 

= tsK-™) 

) (with 1 ++ ), 

0(3900/3885) 

= 7= [DD* + D*D^ 

| (with l 4 ), 

0(4020/4025) 

= D*D (with l" 1 

or 2 ++ ) , 

Y (4140) 

= D*D* S (with 0 ++ ) 
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0(10610) 

= 

(with l 4 ), 

0(10650) 

= B*B* (with 1+-) 



There are more than one possibilities in assigning the 0(3900/3885) and 0(4020/4025), we 
can assign them either in the scenario of tetraquark states or molecular states. The 0(3900), 
0(4200) and Z( 4430) have the same quark constituents ccud or ccdu and analogous decays, see 
Eq.(l). Neither scenario alone can accommodate the 0(3900), 0(4200) and Z( 4430) consistently, 
however, there are rooms to assign those Z c mesons consistently in the two scenarios together, 
see section 4. In Ref. [19], we observe that if we determine the energy scales of the QCD spectral 
densities with the same parameter M c , the 8 g) 8 type molecule-like states have larger masses than 
the corresponding lg)l type molecular states. We obtain the masses M z r 390 o) = 3.891^09 MeV and 
M z c ( 808) = 4.10lg ( j ) o MeV with the QCD sum rules m■ The upper bound of the predicted mass 
M Zc ( 8®8) = 4.10lonoMeV is consistent with the experimental value M Zc ( 4200) = 4196l2gii3 MeV 
PQ. Now, we assign the 0(4200) to be the 8 g) 8 type molecule-like state tentatively, 


Z c ( 4200) 


1 

71 


VV +V*V (with 1+') 


(7) 


study the mass and decay width with the QCD sum rules in details, and fit the effective mass M c 
for the 8 g) 8 type molecule-like states, where the T> and V* have the same quark constituents as 
the D and D* respectively, but they are in the color 8 representation. 

The article is arranged as follows: we derive the QCD sum rules for the mass of the 8 g) 8 type 
axial-vector molecule-like state 0(4200) in section 2; we derive the QCD sum rules for the width 
of the 8 g) 8 type axial-vector molecule-like state 0(4200) in section 3; we discuss the possible 
assignments of the 0(3900), 0(4200) and Z( 4430) as the 3g)3 type axial-vector tetraquark states 
in section 4; section 5 is reserved for our conclusion. 


2 The mass of the 8 08 type axial-vector molecule-like state 


In the following, we write down the two-point correlation function 11^(p) in the QCD sum rules, 


IWp) 


i J d A xe ip ' x (0\T { J M (x) j/(0)} |0), 

u(x)ij 5 \ a c(x)c(x)^ ll X a d(x) + u(x)'y /1 X a c(x)c(x)i^ 5 X a d(x) 


( 8 ) 

(9) 


where the A a is the Gell-Mann matrix in the color space. We construct the 8<8>8 type current J M (x) 
(see Ref.[21] 1 to study the molecule-like state 0.(4200). Under charge conjugation transform C , 
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the current Jp(x) has the property, 

CJ ll (x)C ~ 1 = -Jp(x) | u ++d ■ 


( 10 ) 


We insert a complete set of intermediate hadronic states with the same quantum numbers as the 
current operator J p (x) into the correlation function H pl y(p) to obtain the hadronic representation 
[22l 123] , and isolate the ground state contribution, 


n u, v (p) 


\2 

A Z C ( 4200) 

M 2 — T ) 2 

Z c { 4200) V 



H-, 


(ii) 


where the pole residue Xz c ( 4200 ) is defined by 


(0| J m (0)|Z c (4200)) — A2„(42oo) e p > 


( 12 ) 


the is the polarization vector of the axial-vector meson Z c (4200). 

We carry out the operator product expansion to the vacuum condensates up to dimension-10, 
and obtain the QCD spectral density through dispersion relation, then take the quark-hadron 
duality and perform Borel transform with respect to the variable P 2 = —p 2 to obtain the following 
QCD sum rule, 


a L( 42 oo) ex P ^ - ° 0) j = ^ 2 dsp{s) exp (-^) . (13) 

One can consult Ref.[l9] for the explicit expression of the QCD spectral density p(s). We differen¬ 
tiate Eq.(13) with respect to r = then eliminate the pole residue Xz c ( 4200 ) 1° obtain the QCD 
sum rule for the mass, 


M 2 

1V1 Z C ( 4200) 


SZj ds (~^)p( s ) e ~ 

JC; dsp(s)e- TS 


(14) 


The current Jp{x) has non-vanishing couplings with the scattering states DD*, J/ipir, rj c p , etc 
[24] . In the following, we study the contributions of the intermediate meson-loops to the correlation 
function II^,,^), 


n M „(p) = - 


A i 


p 2 — M| 
p 2 — m| 


■9nv(p) — 9^{p)'P‘DD-{p)g aP {p)gpu{p) Xz ‘ 


p 2 — M§ 


p 2 - M 2 


9 pa (P)£ J/ilm {P)g af) (P)df (p) 


^Z c 


P 2 - M 2 


XZc 9n<x{p)'^r lc p{p)g a0 {p)9p V {p) XZ ‘ 


p 2 — Mj- 


p 2 - M 2 


X Za 


p 2 - M z c - E dd*(p ) - S J/V)7r (p) - T, VcP (p) + • • ■ 


SV(p) + ' 


(15) 


where 


Sdd* (?) = i J 


(2tt)‘ 

^J/' 07 r(p) = i J 


(2tt)‘ 

s Vcp(p ) = i J 


(2 ^ 


C 2 

U Z 0 DD* 


4 \„2 


T 2 1 ’ 


r>2 

^ZcJ/^-K 




r 2 

^ZcT/cP 


[q 2 - M 2 c ] [(p - q) 2 - M 2 \ ’ 


( 16 ) 


4 



















9 ^{p) = + the G Zc dd *, Gz c j/ipm Gz c rj c p are hadronic coupling constants, the X Zc and 

M Zc are bare quantities to absorb the divergences in the self-energies T,od*(p), ^j/^-k{p)j ^i?cp(p)i 
etc. The renormalized self-energies contribute a finite imaginary part to modify the dispersion 
relation, 


n p„(p) 




p 2 — + iy/^Tip*) 


9pv{p) 


(17) 


The present work (see section 3) indicates that T(Z+ (4200) —> J/ijm + ) = 24.6 MeV, T(Z+ (4200) —> 
r] c p + ) = 309.1 MeV and T(Z+(4200) —> (DD*/D*D) + ) = 0 due to the special structure or inter¬ 
polating current of the ^(4200), the width originates dominantly from the decays to the rj c p. We 
can take into account the finite width effect by the following simple replacement of the hadronic 
spectral density, 


( s — 


\2 1 _ M Zc TzA s ) _ 

Zc tt(s- M 2 Za ) 2 + M 2 z V 2 Zc {s) ’ 


(18) 


where 


rz„(s) 


Mzc_ I s - (M Ve + Mtf~ 
s Y M| e - {M Vc + M p ) 2 ' 


(19) 


Then the phenomenological sides of the QCD sum rules in Eqs. (13-14) undergo the following 
changes, 


A| c exp 



Af c M§ c exp 



^z, 


/•So 

J(Mr, 


ds— 


M Zc T Zc (s ) 


+M„) 2 7T (s - M| ) 2 + M§ r| (a) 


exp 


(-£)■ 


0.76(0.79) A| exp - 


Ml 

rp2 


A I, 


ds s -- 


M Zc T Zc {s ) 


(M nc +M p ) 2 7T (s MfJ 2 -f M^T 2 Zc (s) 

9 9 / Af 2 

0.74(0.77) A| c M| c exp(-^f 


exp 


(-?) 


( 20 ) 


( 21 ) 


where we have used the central values of the input parameters, and T^ c = 370 MeV, the values 
0.79 and 0.77 in the bracket originate from the approximation T(Z+) « T(Z+( 4200) — > r] c p + ). So 
we can absorb the numerical factors 0.76 and 0.74 into the pole residue A z c safely, the intermediate 
meson-loops cannot affect the mass M z j 4200 ) significantly, the zero width approximation in the 
hadronic spectral density works. 

We take the standard values of the vacuum condensates (qq) = —(0.24±0.01 GeV) 3 , ( qg s aGq) = 
TO o {qq)i ml = (0.8 ± 0.1) GeV 2 , ( aaGG ) = (0.33 GeV) 4 at the energy scale p = 1 GeV [22) 1221 [25] . 
and choose the MS mass m c {m c ) = (1.275 ± 0.025) GeV from the Particle Data Group [24] . 
Furthermore, we take into account the energy-scale dependence of the input parameters, 


(w) (p) 

(qg s aGq)(p) 

m c (p) 

ot s {p) 


(m)(Q) 


a s (Q) 


( qg s aGq){Q ) 
m c (m c ) 

1 - 


_a s (p,)_ 

a s (Q) 


a s {p)_ 

12 

&s(p) 


1 

b 0 t 


a s (m c ) 

bi log t 6 2 (log 2 1 - log t - 1) + b 0 b 2 


bo t 


bit 2 


( 22 ) 
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T 2 (GeV 2 ) 


Figure 1: The mass of the Z c { 4200) with variations of the Borel parameters T 2 and energy scales 
p,, where the horizontal line denotes the experimental value. 


33 — 2 rif 


7 153—19 rif 

= 24it' 2 ' 


9X^7— 5033 „ | 325 2 

b 2 = £857— , ^Uf+ 27 n f , A = 2 i3 Me V, 296 MeV 


where t = log ^, b 0 = 127r 

and 339 MeV for the flavors rif = 5, 4 and 3, respectively [24]. We tentatively take the continuum 
threshold parameter as ^/so = (4.7± 0.1) GeV 2 , i.e. ^/so = M Zc ( 4200 ) + (0-4 — 0.6) GeV, and search 
for the optimal Borel parameter to satisfy the two criteria (pole dominance and convergence of the 
operator product expansion) of the QCD sum rules. 

In Fig.l, we plot the mass of the Z c (4200) with variations of the Borel parameters T 2 and energy 
scales /r for the threshold parameter y4so = 4.7 GeV. From the figure, we can see that the masses 
decrease monotonously with increase of the energy scales. We can reproduce the experimental 
value M Zc ( 42oo) = 4196±|S±i| MeV [Tj approximately at the energy scale /i = (1.3 — 1.5) GeV. In 
this article, we take the energy scale /.4 = 1.4 GeV. In the Borel window T 2 = (3.0 — 3.4) GeV 2 , 
the pole contribution is about (43 — 64)%, it is reliable to extract the ground state mass. 

In Fig.2, we plot the contributions of different terms in the operator product expansion with 
variations of the Borel parameters T 2 for the threshold parameter ^/so = 4.7 GeV and energy scale 
/i = 1.4 GeV. In the Borel window T 2 = (3.0 — 3.4) GeV 2 , the D 3 D 0 ss |D 5 | D 6 |D 8 |, and 
the D 4 , D 7 and D w play a less important role, where the Di with i = 0, 3, 4, 5, 6 , 7, 8 , 10 denote 
the contributions of the vacuum condensates of dimensions D = i, and the total contributions are 
normalized to be 1. The operator product expansion is well convergent. 

We take into account all uncertainties of the input parameters, and obtain the values of the 
mass and pole residue of the Z c ( 4200), which are shown explicitly in Figs.3-4, 


M Zci 4200 ) = 4.19 ±0.08 GeV, 

\ Zc{ 4200 ) = (5.25 ± 0.71) x 10 - 2 GeV 5 . (23) 

The predicted mass M Zc ^ A 2 00) = 4.19±0.08 GeV is consistent with the experimental value M z c (42oo) = 
4196 I 29-13 MeV within uncertainties [Tj. The QCD sum rules favor assigning the Z c (4200) to be 
the 8®8 type UD* +V*V molecule-like state. Now we can obtain the parameter M c = 1.98 GeV for 
the 8® 8 type molecule-like states according to the energy scale formula /i = J M\/y/z ~ (2M C ) 2 . 
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T 2 (GeV 2 ) 


Figure 2: The contributions of different terms in the operator product expansion with variations of 
the Borel parameter T 2 , where the 0, 3, 4, 5, 6, 7, 8 and 10 denotes the dimensions of the vacuum 
condensates. 



T 2 (GeV 2 ) 


Figure 3: The mass of the Z c ( 4200) with variations of the Borel parameter T 2 , where the horizontal 
line denotes the experimental value. 
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T 2 (GeV 2 ) 


Figure 4: The pole residue of the Z c { 4200) with variations of the Borel parameter T 2 . 

3 The width of the 8 <g>8 type axial-vector molecule-like state 

We can study the strong decays 4200) —> i] c p ± and ( DD *) ± (or (D*D) ± ) with the 

following three-point correlation functions !!„ (p, q), If cL(p, (?) and H^(p,q), respectively, 


n i,(p,q) = i 2 1 d 4 xdW px e iqy (0\T{j^(x)J^(y)M0)}\0), (24) 

n^(p,«) = i 2 J d 4 xd 4 ye^ x e iqy {0\T{J^(x)J^y)M0)}\0), (25) 

n l,(p,q) = i 2 J d 4 xd 4 ye^e iq y(0\T{j°\x)J 5 D (y)J^(0)}\0), (26) 

where the currents 

= c(af )7 a c(x), 

J${y) = u{y)ilsd(y), (27) 

Jn c ( x ) = c(x)ij 5 c(x), 

Ja(y) = u(yhad(y), (28) 

Ja ( X ) = c( X )la.d{x) , 

J°{y) = M(y)i7 5 c(y), (29) 


interpolate the mesons J/if), tt, r] c , p, D* and D, respectively. At the leading order 0(a s ), 
II^ (p, q) = 0 at the QCD side according to the structure of the SU( 3) color group. 

We insert a complete set of intermediate hadronic states with the same quantum numbers as 
the current operators into the three-point correlation functions II* (p, q), n 2 ^^, q) and isolate the 











ground state contributions to obtain the following results, 


^ap(P^ g) 


f. 


r M nfj/4’ M J/ip^Z c G Zc J/^Tr _ 1 _ ( PaPp\ 

m u + m d (M| c - p' 2 ){M 2 ^ - p 2 ){M 2 - q 2 ) \ 9al3 p 2 ) 

_/ + ^ + - 
^ rjl ^ 


= It 


p p' 2 y 


m u + m d 


1 

+ (^i c -p' 2 


1 

(M 2 c -p' 2 )(M 2 ^~p 2 )(MZ- q 2 ) 


1 r ,,pz c *(p\t,p' 2 ) 

P' 2 )(M} H ,-P 2 ) J s 0 t-g 2 


1 

+ (Mf 


z . (gap H-) + • ' 

J3 Z I 


1 Z 100 ,+ Pz c j/4>{t,q 2 ,p' 2 ) 

r | c — p ,2 )( M 2 — q 2 ) t-p 2 

I fnM 2 fj/^Mj/^Xz c Gz c j/4j-K _1_ 

\ m. u + m d (M 2 c - p ,2 )(Mj /ilj - p 2 ){M 2 - g 2 ) 

4_^_+__+ ... 1 ( a - 

+ (M 2 c -p' 2 )(M 2 ^ -p 2 ) + (M 2 c — p ,2 )(M 2 ~q 2 ) + f (9a " 


n « M (fhg) 


/^cZVf 2 fpM p Xz c Gz c ri c p 


2 ?n c 

y+ ^ + .. 

5m p/2 I + 


f Vc MlfpMpX Zc G z , 


(Mz c ~ P' 2 ){M 2 — p 2 )(M 2 — q 2 ) 


Qctfi + 


gqg/A 
g 2 ) 


cVcp 


2m r 


(MZ — P ,2 ){M 2 — p 2 ){M 2 — q 2 ) 


1 

r ^pz cP (P 2 ,t,p' 2 ) 

{M 2 Zc -p l2 ){Ml -p 2 ) 

I UjO - 

t-q 2 

1 

f°° ^Pz c r, c {t,g 2 ,P' 2 ) 

{M‘z c ~ P ,2 ){M 2 — q 2 ) , 

I Ul/- 

/s° 

Vc 

t — p 2 

fr, c M 2 JpMpX Zc Gz cVc p 


1 

2 m c 


p' 2 ){M 2 c — p 2 )(M 2 

Cz cP 

1 

Gz cVc 

(Mz c ~ P l2 )(M 2 c — p 2 ) 

( M l 

— p' 2 ){M 2 - q 2 ) 


(go 


{gap H—) H— 


(31) 


where p' = p + g, the f Vc , / p and f K are the decay constants of the mesons J/V’, Vc, P and 

7r, respectively, the Gz c j/ 4 ,-k and Gz c r] c p are the hadronic coupling constants. In this article, we 
choose the tensor g a p to study the coupling constants Gz c j/ 4 ,-k and Gz c??c p. 
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In the following, we write down the definitions, 


<0|J^(0)|JMp)> 

<0|J£(0)|p(g)> 

(0|^ e (0)|»?c(p)) 

<0|^(0)|7r(g)) 

(J/i/j{p)n(q)\Z c (p')) 

{Vc(p)p{q)\z c (p ')) 


fpMp£ a , 

UK 

2m c ’ 

/^M 2 
m u + m d 

iC(p) ■ ap')G ZcJ/ ^ 

ie*(q) ■ C (p')Gz c r, cP , 


(32) 

(33) 


the £, £ and e are polarization vectors of the J/V’, Z c (4200) and p, respectively. The four un¬ 
known functions Pz c k{p 2 ,t,p' 2 ), Pz c j/ip(t, Q 2 ,p' 2 ), Pz cP (p 2 ,t,p' 2 ) and Pz c r, c {t,q 2 ,p' 2 ) have complex 
dependence on the transitions between the ground states and the high resonances or the continuum 
states. We introduce the notations Cz c k, Cz c j/-4,, Cz cP and Cz c ri c to parameterize the net effects, 


n 

4 

II 

r dt Pz c n(P 2 ,t,p' 2 ) 

Js o t-q 2 


GZaJ/lp ~ 

f°° ^ Pz c j/ip {t,q 2 , p ' 2 ) 
f - p2 


n 

II 

Js° t - q 2 


Cz cVc = 

JsO t - P 

(34) 


In numerical calculations, we smear the dependencies of the Cz c7r , Cz c j/ip, Gz cP and Cz c rt c on the 
variables p 2 , p' 2 , q 2 , take the Cz c7r , Cz c j /^, Cz cP and Cz cVc as free parameters, and choose the 
suitable values to eliminate the contaminations to obtain the stable sum rules with the variations 
of the Borel parameters [26j. 

We carry out the operator product expansion up to the vacuum condensates of dimension 5 
and neglect the tiny contribution of the gluon condensate, one can see Fig.2 as an example. In the 
QCD sum rules involving the tetraquark states, if there exist contributions from the perturbative 
terms or quark condensate terms, or mixed condensate terms, then the gluon condensates play a 
minor important role. The leading-order perturbative terms, the quark condensates, the mixed 
condensates and the gluon condensates are the vacuum expectations of the operators of the order 
G(a°), O{op s ), 0(«y 2 ) and C?(a s ), respectively, so the gluon condensates can be neglected approx¬ 
imately. Furthermore, the double Borel transformed QCD sum rules converge much faster than 
the single Borel transformed QCD sum rules at the operator product expansion side. We obtain 
the QCD spectral densities through dispersion relation, take the quark-hadron duality below the 
continuum thresholds, then set p' 2 = p 2 and take the double Borel transforms with respect to the 
variable P 2 = —p 2 and Q 2 = —q 2 respectively to obtain the following QCD sum rules, 


fTT M nfj/ip^J/ip^Z c G Zc J/^T l 


m u + md 


M k - M jh, 


exp - 


M W 

T 2 


Ml 

- ex P( —ff 


exp 


Ml 


Ml 


+Cz a j/i, exp ^ 

m c (qg s aGq) 


12V27T 4 


ds 


duu (s + 2m 2 ) \ 1 ^ 


4 ml 


ex P( -jq ~fa 1 ~ 


6\Z2tt 2 


' 4 


ds\ 1 — 


4m? 


ex P I - fa 



(35) 
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fvc M Vc fp M p^z c Gz c 


VcP 


2 m c 

+Czcr/c ex P 
m c (qg s aGq) 


+ 


18 a/ 2 ti 


n 


U ml 


M 2 Zc -Ml 
Mf\__ 

n 


exp 


12yph r 4 Jim? 


K 

T ? , 

Zc ds 


exp 


Ml 

T? 


exp 


4 777,2 

duus\/l -- exp 


ds\ 1 — 


4m? 


ex P 1 




(36) 


where the s° Zc , s 6 and s 6 are the continuum threshold parameters for the Z c (4200), 7r and p, 
respectively. 

The hadronic parameters are taken as M= 0.13957 GeV, M p = 0.77526 GeV, Mj= 
3.0969GeV, M Vc = 2.9836GeV [24], = 0.130GeV, f p = 0.215 GeV [25], fj^ = 0.418GeV, 

f ric = 0.387GeV [27], = 0.85GeV, = 1.3GeV [25], = 4.7GeV, A Zc = 5.25 x 

10~ 2 GeV 5 , T 2 = (3.0 — 3.4) GeV 2 (present work), T 2 = (0.8— 1.2) GeV 2 [25] . /,Mj/(m„ + 
rrid) = —2 (qq)/fn from the Gell-Mann-Oakes-Renner relation. The unknown parameters are cho¬ 
sen as Cz c j/ij) = 0.01 GeV 6 and Cz cVc = 0.09 GeV 6 to obtain platforms in the Borel windows 
T 2 = (3.0 — 3.4) GeV 2 . The parameters at the QCD side are chosen as the same in the two-point 
QCD sum rules for the Z c (4200). Then it is easy to obtain the values of the hadronic coupling 
constants, 


Gz c j/if>v = 3.34 ±0.07 ±0.25 GeV, 

Gz cVcP = 11.31 ± 0.07± 1.06GeV, (37) 

where the uncertainties originate from the Borel parameters T 2 and T 2 , respectively, see Figs.5-6. 
As the largest uncertainties originate from the Borel parameter T 2 , we neglect the uncertainties of 
the parameters other than the Borel parameters. The uncertainties of the Gz c j/^n and Gz c r/ c p lead 
to the uncertainties 5T(Z+(4200) —> JM7r + )/T(Z+(4200) -» J/itm + ) = 2 SGz j/^-k/Gz j/p-k ~ 
8 % and flW(4200) -> „ WtA 4 \j) = ^z.JlGzl * 9%. 

The central values of the decay widths are 


T(Z+(4200) -± J/ip 7T+) 


T(Z+(4200) -> r, c p+) 


V (Mz c , , Ml) 2 0 P (M Zc , Mj , Ml) 

2A-kM 2 Zc z ‘ j I^ + Mj /q/> 

24.6 MeV, 


P ( Mz c , M rjc , M p ) 2 o . P {Mzc i M Vc , M p ) 
24tt M 2 Zc z ^ p M 2 

309.1 MeV, 


(38) 


where b, c) = 


2a . If we saturate the width of the Z c (4200) with the strong 

decays to J/ipn and p c p, then T Zc ( 42 oo) ~ 334 MeV, which is consistent with the experimental 
value r Zc ( 4200 ) = 3704 oti 32 MeV from the Belle collaboration [T], the present calculations support 
assigning the Z c (4200) to be the 8 ® 8 type axial-vector nrolecule-like state. Due to the special 
structure of the 2), (4200), the decays to the final states DD* and D*D can only take place through 
re-scattering mechanism Z c (4200) —> J/ij)n, r] c p —> DD*, D*D , the decay widths T(Z C (4200) —> 
DD*, D*D) are expected to be small. 
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Figure 5: The coupling constants Gz c j/^n {A) and Gz^p (B) with variations of the Borel pa¬ 
rameters Tf. 
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Figure 6: The coupling constants G Zc j/^-k {A) and Gz cVcP (B) with variations of the Borel pa¬ 
rameters T% ■ 


12 





















4 The masses of the 303 type axial-vector tetraquark states 

In the following, we write down the two-point correlation function II^^p) in the QCD sum rules, 


rV(p) = iJd^xe^iOlT^ixHi 0)}|0), (39) 

f ijk Amn 

V^{x) = — —J= — {u J (x)C'y 5 c k (x)d m (x)'Y IJ ,Cc n (x) - u J (x)Cj li c k (x)d m (x)'y 5 Cc n (x)} ,(40) 

the *, j, fc, m, n are color indexes, the C is the charge conjugation matrix. We choose the current 
r)^{x ) to interpolate the 3 g) 3 type axial-vector tetraquark states. Under charge conjugation 
transform C', the current 70 a;) has the property, 

Ci] tl (x)C~ 1 = -70a;) Iu-H-d ■ (41) 


We carry out the operator product expansion up to the vacuum condensates of dimension 10 and 
obtain the correlation function at the QCD side, 


10(0 



ds 


P(s) 

9 

S — p z 



PliVv\ 

p 2 J 


Tf / 2 \P^Pv 

n 0 (p )-0- 


(42) 


the expression of the QCD spectral density p(s) is shown explicitly in Ref.|17]. the component 
IIo(p 2 ) is irrelevant in the present analysis. 

In case I, the 0(3900) and 04430) are the ground state and the first radial excited state 
of the 3 (8) 3 type axial-vector tetraquark states, respectively, the 0(4200) is not the 3 (g> 3 type 
axial-vector tetraquark state, then the current couples potentially to the 0(3900) and 04430), 
(O|t 0O)|0(39OO)/0443O)) = Az c ( 3900 )/.z( 4430 ) e Ao where the e M are the polarization vectors of the 
0(3900) and 04430). Now we retain the ground state and the first radial excited state and write 
down the QCD sum rule m , 


\2 ( M Z c (3900)\ „( M Z{ 4430)^ 

0(3900) ex P I T 2 ) ^ z ( 443 °) exp 1 T 2 ) 



We differentiate Eq.(43) with respect to r = 0 and obtain three additional QCD sum rules, 


A 2 1\/T 2 I zj c y*j&\j\jj i » ^ 1\/T 'Z I zj y-t^ 

Z c (3900) M Z C (3900) ex P I ^2 I + * Z(4430) M Z(4430) ex P I ^2 


M 2 

111 Z c ( 3900) 


M 2 ' 

iW Z(4430) 


/ 4m^ 


ds sp(s) exp (-^ 5 ) , 


(44) 


M 2 

1V1 Z C ( 3900) 


\ Z ]\ I zj c yu 

A Z C (390000 (3900) ex P 1 ^2 

ds s 2 p(s) exp (-^j) 


rs 0 

/ 4m^ 


M 2 s 
1V1 Z{ 4430) 


+ ^Z(4430)Mz (4430) ex P I ^2 


(45) 


A| c (3900) M Zc (3900) ex P 
so 


M 2 

1V1 Z C ( 3900) 


0*,V(»)exp(-^) . 


M 2 

JU Z(4430) 


4 ”") 1 \ 2 1/6 

Tpi I + A Z(4430) JW Z(4430) ex P 1 7^ 


(46) 
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We solve the coupled equations consistently and obtain the values of the masses of the £ c (3900) 
and £(4430), 

Mz c (. 3900 ) = 3 . 9 lio!i 7 GeV, Experimental value 3899.0 ± 3.6 ± 4.9 MeV [5] , 

Mz( 4430 ) = 4.5llo;Jg GeV, Experimental value 4475 ± 7 +25 MeV [4] , (47) 

which favors assigning £ c (3900) and £(4430) to be the ground state and the first radial excited 
state of the 3 <g> 3 type axial-vector tetraquark states, respectively. For the technical details, one 
can consult Ref. [12] . We can assign the £ c (4200) to be the 8 ® 8 type axial-vector molecule-like 
state, or it is possible to assign the £ c (4200) to be the 8 g> 8 type axial-vector molecule-like state. 

In case II, the £ c (4200) is the ground state of the 3 ® 3 type axial-vector tetraquark state, the 
£ c (3900) and £(4430) are not the 3g)3 type axial-vector tetraquark states, then the current couples 
potentially to the £ c (4200), (0|?7^(0)|£ c (4200)) = ^z c (42oo) £ v’ where the £ n is the polarization 
vector of the £ c (4200). Now we write down the QCD sum rule, 

^i c ( 4200 ) exp ^- Zc ^f° 0) j ^dspis) exp . (48) 

We differentiate Eq.(48) with respect to r = j %■, then eliminate the pole residue A^ c ( 42 oo) to obtain 
the QCD sum rules for the mass, 


A-/ 2 

IU Z C (4200) 


IZ* da (- 


) P(s)t 


JZ? dsp{s)e~ 


(49) 


In Fig.7, we plot the mass of the £ c (4200) with variations of the Borel parameters T 2 and energy 
scales p for the threshold parameter = 4.7 GeV. From the figure, we can see that the mass 
decreases monotonously with increase of the energy scales, the experimental value Mz c ( 4200 ) = 
4196±^i^MeV [T] can be reproduced approximately at the energy scales p = (1.1 — 1.4) GeV. If 
the £ c (4200) is the ground state of the 3 ® 3 type axial-vector tetraquark state, we have to assign 
the £ c (3900) to be the 1 ® 1 type molecular state [l9|, 


£ c (3900/3885) 


-L(bd- + d>d), 


(50) 


however, it is odd to assign the £(4430) to be the excited 1 (g) 1 type molecular state. 

In case III, the £(4430) is the ground state of the 3 ® 3 type axial-vector tetraquark state, 
the £ c (3900) and £ c (4200) are not the 3 (g) 3 type axial-vector tetraquark states, then the current 
couples potentially to the £(4430), (0|?y /i (0)|£(4430)) = Xz( 4430 ) £ t^ where the is the polarization 
vector of the £(4430). Now we write down the QCD sum rule, 


A l( 4430 ) ex P ^ - Z y 2 430) j ^ ^dsp(s) exp(-^j) . (51) 

We differentiate Eq.(51) with respect to r = eliminate the pole residue \z( 4430 ) to obtain the 

QCD sum rules for the mass, 


M 2 

m Z( 4430 ) 


IZi ds Z-t)p( s ) e ~ 

IZ’ dsp(s)e~ TS 


(52) 


In Fig. 8 , we plot the mass of the £(4430) with variations of the Borel parameters T 2 and energy 
scales p for the threshold parameter ^/so = 5.0 GeV. From the figure, we can see that the experi¬ 
mental value M Z ( 443 o) = 4475 ± 7 ^25 MeV [3] can be reproduced approximately at the energy scale 
p = 1 GeV [33]. At the energy scales p > 1 GeV, the predicted mass M Z ( 443 o) is much smaller than 
4475 MeV, assigning the £(4430) to be the ground state of the 3 (g) 3 type axial-vector tetraquark 
state is not favored. 
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Figure 7: The mass of the Z c { 4200) with variations of the Borel parameters T 2 and energy scales 
where the horizontal line denotes the experimental value. 



T 2 (GeV 2 ) 


Figure 8: The mass of the Z( 4430) with variations of the Borel parameters T 2 and energy scales 
where the horizontal line denotes the experimental value. 
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5 Conclusion 


In this article, we assume the ^(4200) as the 8 0 8 type axial-vector molecule-like state, and 
construct the 8 0 8 type axial-vector current to study its mass and width with the QCD sum rules. 
The numerical result supports assigning the Z c ( 4200) to be the 80 8 type molecule-like state with 
J PC = l" 1 . Furthermore, we discuss the possible assignments of the Z c ( 3900), Z c ( 4200) and 
Z( 4430) to be the 3 0 3 type tetraquark states with J pc = l -1 . 
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